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We simulate the optical and electrical responses in gallium-doped graphene. Using density func¬ 
tional theory with a local density approximation, we simlutate the electronic band structure and 
show the effects of impurity doping (0-3.91%) in graphene on the electron density, refractive index, 
optical conductivity, and extinction coefficient for each doping percentages. Here, gallium atoms 
are placed randomly (using a 5-point average) throughout a 128-atom sheet of graphene. These 
calculations demonstrate the effects of hole doping due to direct atomic substitution, where it is 
found that a disruption in the electronic structure and electron density for small doping levels is 
due to impurity scattering of the electrons. However, the system continues to produce metallic or 
semi-metallic behavior with increasing doping levels. These calculations are compared to a purely 
theoretical 100% Ga sheet for comparison of conductivity. Furthermore, we examine the change in 
the electronic band structure, where the introduction of gallium electronic bands produces a shift 
in the electron bands and dissolves the characteristic Dirac cone within graphene, which leads to 
better electron mobility. 


I. INTRODUCTION 


Understanding the complex properties of materials is 
a fundamental goal for the advancement of technologies 
and other industries in the world today. Specifically, the 
examination and prediction of electronic properties and 
phase transitions in insulators and metals has gained 
much attention, typically for the investigations of ex¬ 
otic topological states and/or large symmetry-breaking 
regimes^. Topological materials in particular have been 
gaining a lot of attention due to their many surface prop¬ 
erties. Specifically, the study of the electronic properties 
of graphene and its interactions with substrates and other 
materials has show to exhibit many desired phenomena 
for technological devices^. 

Graphene is a sp^-bonded 2D sheet of carbon atoms 
arranged in a honeycomb pattern (shown in Fig. . 
The bonding and arrangement of atoms allows graphene 
to have a very large tensile strength^ and interesting 
thermal properties'*^, which when considering its size 
and weight, is stronger than steel^. Furthermore, the 
TT-bonding throughout the lattice provides a high elec¬ 
tron mobilit} ^^^^ , where the high durability, strength, 
and conductivity make graphene a suitable choice for use 
in nanodevices and electronics. However, even though 
graphene has a high tensile strength, the ability for the 
material to be a great component for electronics is still 
unclear. 

Currently, Graphene is being investigated for use in 
electronics through surface deposition on metals and 
semiconductors, where the interactions of graphene and 
other atoms is hoped to produce enhancements to the 
overall electronic properties'^. This enhancement is due 



FIG. 1: The structure (left) and calculated energy density 
(right) for a pure graphene sheet (top) and theoretical gallium 
sheet (bottom). 


to the presence of surface or topological states^, which 
are critical to the understanding of the electron mobility 
in graphene, and has led to many researchers working to 
find ways for improving the electronic structure and over¬ 
all electronic properties of graphene. While topological 
interactions with substrates and other materials plays an 
important role these advancements, the ability to change 
and control the electronic properties of graphene through 
direct substitution could provide a unique tunability to 
the mechanisms that drive technology further. 
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FIG. 2: The calculated electron density as a function of spatial position in graphene with various concentrations of gallium 
(from left to right - 0.78, 1.56, 2.34, 3.13, and 3.90% Ga doping). The crystal structure is given directly below the density 
plots. 


Recent experiments have demonstrated that it is possi¬ 
ble to remove individual carbon atoms from the graphene 
latticei^. Furthermore, various experimental and theo¬ 
retical papers have also shown that carbon atoms can 
be replaced with other atoms (i.e. nitrogen, cobalt, and 
iron) in specifically designed spotsi^^— . This opens a 
whole realm of possibilities for the tuning of the elec¬ 
tronic properties of 2D materials through a manipulation 
of the electronic and topological states. 

Topological states in graphene are produced through a 
unique conductivity point in its band structure known as 
the Dirac point, which typically allows electrons to easily 
move between the conduction and valence bands through 
a band structure Dirac cone at the Fermi level. The in¬ 
troduction of substituted or doped atoms will produce 
a breaking of the Dirac cone. However, the presence of 
these atoms will produce an impurity resonance at the 
Fermi level, which should provide an increase in the elec¬ 
tron mobility and overall conductance. In many case 
thus far, researchers have examined electron favorable or 
magnetic atoms. Therefore, we propose the idea that 
impurity states may be hindered through hole doping or 
p-type substitution in graphene using gallium atoms. 

To examine the effects of topological states produced 
by p-substituted atoms, we examine the presence of 
directly-substituted Gallium atoms in graphene at low 
doping levels of 0 - 3.91%. Using density functional the¬ 
ory in an local density approximation (LDA), we deter¬ 
mine the electronic structure and optical properties for 
the various stages of Ga doping. Overall, we find that 
low doping of graphene with Ga will actually produce 
a decrease in the overall conductivity, which is a trend 
observed in the index of refraction and extinction co¬ 


efficient as well. To reconcile these observations with 
known production of impurity resonances in Dirac ma¬ 
terials, we examine the electron density and compare 
calculations to a theoretical 2D structure of pure Ga 
atoms. Here, we determine that the presence of a small 
amount p-type dopant produces impurity scattering in 
the graphene sheet, which is the likely cause for the de¬ 
crease in overall electron mobility. 

II. COMPUTATIONAL METHODOLOGY 

We ran a series of density functional calculations using 
Atomistix Toolkit by QuantumwiseiSr— . In this study, 
constructed a 128 atom supercell of graphene and sub¬ 
sequently replace random carbon atoms with Ga using 
doping levels of 0, 0.78, 1.56, 2.34, 3.125, 3.91, and 100%. 
This is illustrated in Fig. [21 To simulate unbiased place¬ 
ment, Ga atom positions were determined by a random 
number generator, and for systems with 2-5 impurity 
atoms, five different configurations were used to find the 
average quantities. 

Using a local density approximation (LDA), a geom¬ 
etry optimization and band structure calculations were 
performed using a k-point sampling of 5x5x1 with pe¬ 
riodic boundary conditions for 23 different configura¬ 
tions. Furthermore, we determined the optical spectrum, 
and electron density and calculated the conductivity, re¬ 
fractive index, and extinction coefficient. The doped 
graphene structure was kept in the planar configuration 
to allow for us to examine only the electronic effects of 
substitution without energy parameters provided by lat¬ 
tice distortions. Further studies will analyze how the 
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FIG. 3: The simulated electronic band structure for graphene sheets with 0%(a), 0.78%(b), 1.56%(c), 2.34%(d), 3.13%(e) 
3.91%(f), and 100%(g) substitution of gallium. 
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FIG. 4: Simulated dielectric response as a function of energy 
for a pure graphene sheet and a theoretical gallium sheet. 
The bottom panel shows the calculated refractive index and 
extinction coefficient as functions of energy, where the con¬ 
ductance is placed in the inset. 


energy bands, conductivity, and optical properties are 
affected in the presence of a lattice distortion. The ge¬ 
ometry optimization did produce a planar distortion due 
to the presence of the larger Ga atoms. 

From the states in the electronic band structure and 
electron density^!, the Kubo-Greenwood formula can be 
used to calculate the susceptibility tensor from state m 
to n is 


Xii(w) 


m'^eoVuj'^ 


E 


f{Em) - fjEr,) 
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< n\ja\m >< rn\jp\n >, 


( 1 ) 

where ja describes the momentum operator between ma¬ 
trix element between state n and m, V is the volume, rj is 
a finite broadening, and f{Em) = l/(exp(i?/fc_Br)-|-l)2^. 
From the susceptibility tensor—, we determine the real 
and imaginary parts of the dielectric tensor 


er(a;) = (l-hx(w)) (2) 

and the optical conductivity 

cr(a;) = (3) 

From the dielectric tensor—, the refractive index and ex¬ 


tinction coefficient can be calculated as 


+ ^2 + 


(4) 


and the extinction coefficient 


\/^i + ^2 ~ 


(5) 


Using these quantities, it is possible to evaluate the 
electron mobility of supercell. 


III. RESULTS AND DISCUSSION 

Figure [3] shows the calculated electronic band struc¬ 
ture for each configuration. This illustrates the dramatic 
change in electron mobility as Ga atoms are substituted 
into the graphene lattice. Once one Ga atom is substi¬ 
tuted into graphene (Fig. ISKb)), the characteristic Dirac 
cone shifts in energy and is broken and produces a dis¬ 
tinct bands at the Fermi level. As more atoms are substi¬ 
tuted into the 2D lattice, the impurity bands form at zero 
energy which produces the observed impurity resonances. 

From the calculated band structure, the dielectric ten¬ 
sor can be determined as a function of energy. Figure0Ka) 
shows the determined optical spectra for pure graphene 
and a theoretical analog of Ga atoms. This illustrates the 
expected increase in conductivity, which is characteristic 
of the simulated band structure. 

Figure [S] shows the real and imaginary components 
of the calculated dielectric tensors for of doping cases. 
The dielectric response shows a low energy shift as Ga is 
added to the system. From the calculated optical spec¬ 
tra, we can determine the conductivity, refractive index, 
and extinction coefficient (shown in Fig. |6l). 

From the comparison of the optical spectra properties 
for graphene and low doping cases of Ga, there is a dis¬ 
tinct lose of conductivity and electron mobility, which is 
in contradiction to the change in the optical properties 
from graphene to a pure Ga sheet. The optical properties 
also appear to be in slight contradiction to the electronic 
structure for the low doping Ga cases, where the impurity 
bands seem to provide cross-over and conduction states 
at the Fermi level. 

To understand these contradictions, we examined the 
electron density (shown in Figs. [T] and [2]) correlating 
to each configuration. From the electron density in the 
graphene sheet, it is shown that the electrons are dis¬ 
tributed fairly evenly in between each carbon atom in the 
pure graphene sheet, where blue regions represents no or 
a low electron density and red correlated to high elec¬ 
tron density. However, when the first dopant is added, 
the electron density shifts away from this region, which 
is indicated by the large blue region. This shift indicated 
the movement of the electrons from the gallium atoms. 
This trend continues through the low concentration Ga- 
doped sheets. However, in the theoretical pure gallium 
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FIG. 5: Simulated i-axis dielectric response for the vari¬ 
ous concentrations of gallium (0, 0.78, 1.56, 2.34, 3.13, and 
3.91%). This is the average dielectric response consisting of 
a five-point configuration average. The average dielectric re¬ 
sponse along the y axis is similar. 


sheet, the electrons tend to cluster right on top of each 
gallium atom rather than in between the atoms, which 
produces the large increase in the electronic mobility. 

As we have seen in both the optical spectrum and elec¬ 
tron configurations, there is a general disturbance in the 
system causing a decrease in electron mobility. From the 
electron density, this decrease in conductivity seems to 
be due to impurity scattering at the low doping levels, 
which will restrict the ability to enhance the electronic 
properties through direct substitution. The formation 
of low electron density holes from the Ga atoms produce 
transmission pathway disruptions, and as a gallium atom 
is added to the system, electrons encounters the holes in 
the graphene sheet and is scattered from them. This 
trend continues for all low doping percentages and the 
only change seems to occur at the pure gallium sheet. 
In this theoretical gallium sheet, the optical spectrum 
shows an extremely high conductance due to the disap¬ 
pearance of the impurity scattering from the uniformity 
of the sheet. The low to high conductivity from the dop¬ 
ing cases to the theoretical Ga sheet indicates there must 
be a cross-over point in the intermediate region of doping. 
However, when that region was calculated the 2D sheet 
was found to be unstable in the geometry optimization. 
Therefore, it is only possible to examine the low-doping 
regions to understand the possibility of enhanced elec¬ 
tronic and optical properties. 



FIG. 6: The calculated conductance (a), refractive index (b), 
extinction coefficient (c) as function of energy for the various 
substitutions concentrations of gallium (0, 0.78, 1.56, 2.34, 
3.13, and 3.91%). 


tione, it was found the as Ga dopants are added to the 
system, there is a distinct shift in the optical spectrum 
that indicates a decrease in the conductance. Through 
an analysis of the electron for all configurations, it ap¬ 
pears this is due to impurity scattering at the Ga atomic 
site. This analysis was done in the confines of the planar 
geometry to minimize contributions produced by geomet¬ 
ric distortions. Further work will incorporate the z-axis 
graphene distortion to provide a better understanding of 
the system. 


IV. CONCLUSIONS 

In this study, we have investigated the effect of low 
percentage gallium doping on graphene in order to un¬ 
derstand the possibility of electronic enhancement. Using 
density functional theory and a local density approxima- 
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